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With the latest astronomical data including Cosmic Microwave Background (WMAP three year, 
CBI, ACBAR, VSA), Type la Supernova ("gold sample"). Galaxy Clustering (SDSS 3-D matter 
power, Lyman-a forest and Baryon Acoustic Oscillating (BAO)), we make a global fitting to con- 
strain the mass varying neutrinos. We find that the parameter 5, denoting time evolving of neutrino 
mass, is weakly constrained and the neutrino mass limit today can be relaxed at least by a factor of 
two. Adding data of Oi^2/3 decay of Heidelberg-Moscow experiment to our analysis, we find that 5 
can be constrained tightly and mass varying neutrinos are favored at about 99.7% confidence level. 



. . . The mass varying neutrino which has been studied actively in the literature recently brings forth new theoretical 
challenges and abundant phenomenology. The predictions of the variation of the neutrino masses can be tested in 
the experiments of neutrino oscillation Ij, short gamma ray bursts [2] and extremely high energy cosmic neutrino 
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Cosmologically the mass varying neutrino plays an interesting role in determining the evolution of the universe and 
, . has interesting implications in leptogenesis[H, Supernovaf^, the Cosmic Microwave Background(CMB) and Large 
J3 ; Scale Structure(LSS);7]. 

i Neutrino mass variation can be induced by the interaction between neutrino and scalar field [l,[^,[l3,[lH in models 
C ' of neutrino dark energy^. In the framework of ACDM the variation of neutrino mass could result from the coupling 
' . between neutrino and Ricci scalar in the form of f{R)vv and possibly some other unknown mechanism. In this paper 
we work in a model-independent way to constrain the time evolving of neutrino mass by parameterizing neutrino mass 
^ . in the form of Eq.([T]). For simplicity and to extract the effect of mass varying neutrinos we use the simplest dark 
— ' energy model, the cosmological constant for the analysis in this paper. 
04 We take a global fit to constrain the parameter 5 denoting the mass variation of neutrinos and study its effect on the 

determination of other cosmological parameters. Our result shows that l)with the mass variation, the cosmological 
limit on the neutrino mass is relaxed by roughly a factor of two; 2)if taking the limits on the neutrino mass from 
Heidelberg-Moscow experiment (HM) as a prior [3, 15, [l3|, we find that the variation of neutrino mass is tightly 
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; constrained and the mass varying neutrinos are favored at about 3 a. 
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We assume that three species of neutrinos whose masses are degenerate and parameterize the neutrino mass in form 

of 2 



O' m^(a) = m^o[l + ^(1 - a)] (1) 

' where m^{a) is the sum of neutrino masses which is a function of scale factor a, m^o is the sum of neutrino mass at 
, present, d is the dimcnsionless factor denoting the time-varying effect of neutrino masses. For the background, 

P^ = ^ f q^dqdnefoiq), (2) 



M with — + rnl(a)a'^., = ap* is the comoving momentum. The pressure is 



Thus, 



If q 
Pi^ = / q'^dqdQfoiq) — . (3) 



+ 3H{p, + p,) = -^^a(p, - 3p,). (4) 



^ In Ref. '125 Afshordi et al. concluded the neutrino dark energy model presented by Fardon et al.^ is unstable due to the negative sound 
speed squared. However, their argument is based on the ad hoc model in J)j thus lack the generality. Generically, for the mass varying 
neutrino model originally proposed in 8], there exists large parameter space where the system is stable. For a very recent survey, see 
Ref.lll. 

^ The physical motivation for mass varying neutrino mainly stems from the coupling of neutrino and the dark energy thus we choose 
this parametrization which of similar form of widely-used parametrization for dark energy, namely, w(a) = wq -|- uii(1 — a). Here m^o 
and 5 denote the neutrino mass presently and its derivative with respect to scale factor, say, dm{a) / da. In this sense, we work in 
a model-independent fashion. We stress that our parametric method using Eq.l[TJ is just one possible way to study the variation of 
neutrino mass. With the accumulation of high quality astronomical data, we can investigate the neutrino mass variation in a more 
model-independent way, such as using the technique of Principle Component Analysis{PCA). 
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where the overdot represents the derivative with respect to the conformal time. For the evolution of perturbation, 
working in the synchronous gauge and following the treatment of Ref. we find that the Boltzmann equations 
Eq.(40) in 17| don't relate directly to the mass varying terms dlnm^/da, and therefore the hierarchy equations for 
the massive neutrinos Eq.(56) in 17[ remain unchanged. In our numerical calculation we integrate these hierarchy 
equations rather than use the approximate strategy suggested in [3] to assure a high precision. Using a modified 
code of CosmoMc[3|, we make a global fit to constrain the mass varying neutrinos. Assuming a flat universe, we set 
the most general cosmological parameter space as: 



P = (wb, ujc, Qs, T, /^, 6, Us, log[10^%]) , 



(5) 



where uji, = flhh^ and Uc — flch^ are the physical baryon and cold dark matter densities relative to the critical density, 
Qs characterizes the angular scale of sound horizon, r is the optical depth to the last scattering surface and is 
defined as the amplitude of initial power spectrum, fi, is the dark matter neutrino fraction at present, namely. 



PDM 93.105 eV Qch^ 



(6) 



where pi, and pdm denote the energy density of neutrino and dark matter at present respectively. We sample in 
above an eight dimensional parameter space and fit the theoretical output to the observation using Markov Chain 
Monte Carlo algorithm 24, H,|2^. We take the weak priors as: r < 0.8, 0.5 < rig < 1.5, < < 0.5, a cosmic age 
tophat prior as 10 Gyr< to <20 Gyr. To keep the positivity of neutrino mass we take — 1 < (5 < 10. Furthermore, 
we make use of the Hubble space telescope (HST) measurement of the Hubble parameter Hq = IQQh km s~^Mpc~^ 
(sH by multiplying the likelihood by a Gaussian likelihood function centered around h = 0.72 and with a standard 
deviation cr = 0.08. We impose a weak Gaussian prior on the baryon density Ut,h^ = 0.022 ± 0.002 (Icr) from the Big 
Bang nucleosynthesis [ill. 

For CMB data, we use the three-year WMAP (WMAP-3) Temperature- Temperature (TT) and Temperature- 
PolarizationfTE) power spectrum with the routine for computing the likelihood supplied by the WMAP team'^ 
10, mi, m Hi as well as ACBAR CBI and VSA ^ data. To break the degeneracies of the cosmological 

parameters, we add non-CMB data into our analysis. For supernova type la (SN la) of "Riess gold sample" [33|, we 
have marginalized over the nuisance parameter |35l in the calculation of SN la likelihood. For LSS information, we 
have used the 3D matter power spectrum of SDSSfST*] and 2dFGRS [3^, Lyman-a forest data (Lya) from SDSS [3§| 
and the recent measurement of the baryon acoustic oscillation (BAG) feature in the 2-point correlation function of 
SDSS [401. To be conservative but more robust, we only use the first 14 bins of the SDSS 3D matter power spectrum, 
which are well within the linear regime 41 1. For Lya likelihood, we modify the interpolating code^ to incorporate our 
models. For BAG likelihood, we use the constraint [iol : 



A = L>y(0.35) 



0.35c 



0.469 ±0.017 



(7) 



Dviz) = 



DM{zf 



cz 



H{z) 



1/3 



(8) 



where H[z) is the Hubble parameter, c is the speed of light and Dm{z) is the comoving angular diameter distance 
at a specific redshift z. Moreover, the Heidelberg-Moscow experiment uses the half time of Oz/2/3 decay to con- 
strain the effective Majorana mass and this translates to the constraint on the sum of neutrino masses under some 
assumptions [42|: 



m^a - 1.8 ±0.6 eV (2cr) 



(9) 



Given that the Heidelberg-Moscow experiment is controversial for the time being, we just make a tentative fit choosing 
the HM prior. 

For each regular calculation, we run 6 independent chains comprising of 150,000-300,000 chain elements and spend 
thousands of CPU hours to calculate on a cluster. The average acceptance rate is about 40%. We eliminate the first 



^ Available at http: / /lambda. gsfc.nasa.gov/ product / map / current / 1 

* Available at ^http:/ /www.cita.utoronto.ca'^pmcdonal/LyaF/public.lyafchisq.tar.gz] 
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TABLE I: Mean and la constraints on the cosmological parameters. For the weakly constrained parameters, such as 
TTii/o and S for some data combinations, we quote the 95% upper hmits instead. Upper part of the table is for ACDM 
+ neutrinos with constant mass while in the lower part we free 5 to study the mass varying neutrinos. "ALL" denotes 
WMAP3+ACBAR+VSA+CBI+RIESS+SDSS+2dF+Lya throughout this paper. 



5 = ALL ALL+BAO ALL+BAO+HM 

m^o (eV) < 0.616 < 0.393 0.760l:°o°ill 

0.317 ±0.021 0.280 ±0.015 0.303to.M? 
as 0.832 ± 0.024 0.834 ± 0.024 0.795tg g2^ 
S free ALL ALL±BAO ALL±BAO±HM 

S < 6.661 < 8.562 < -0.713 

m^o (eV) < 1.619 < 0.776 l.dGSto]^^ 

n„ 0.319 ±0.024 0.281 ±0.014 0.298 ±0.015 
as 0.829 ±0.027 0.835lo;o24 0.790lo;o23 
Ax^ 0.198 0.222 9.752 



10% chain elements for "build in", and for the convergence test typically we get the chains satisfy the Gelman and 
Rubin 43] criteria where R— 1 < 0.1. 

We summarize our main results of the mass varying neutrinos in the lower part of Table L For comparison, we also 
study the ACDM model with neutrinos of constant mass. The first discovery is that the neutrino mass limit at present 
epoch can be relaxed dramatically if neutrino mass varies during evolution. Without Heidelberg-Moscow data, we 
find that the neutrino mass limit can be relaxed by a factor of 2.6 for "All" data and 1.8 for "AU+BAO"^. Adding 
Heidelberg- Moscow prior, we find that the mean value of neutrino mass rises up from 0.760 to 1.568 while reducing 
by 9.752. This is expected since the Heidelberg-Moscow prior is in great tension with the cosmological observations. 
For example, the authors of Ref. [i^ argued that given the current cosmological constraint on the (constant) neutrino 
mass the HM prior can be excluded. However, this controversy can be resolved if neutrino mass varies. For our ad hoc 
parametrization ([T]), we see that the "All+BAO+HM" prior can put stringent limit on 5, namely, —1<S< —0.713. 
The best fit value of 6 is about —0.9. It means that, in order to be consistent with all the cosmological observations, 
the neutrino mass must be very small in the past, but has grown recently in order to agree with the Heidelberg-Moscow 
data. 




FIG. 1: One dimensional posterior distribution of neutrino mass today mo and its time variation 5. Solid curves denotes the 
mass varying models (see text) while dashed lines show the models with constant neutrino mass. Different data combinations 
are distinguished by color. Green: All; Blue: A11±BA0; Red: AU±BAO±HM prior. 



Mass varying neutrinos lessen the tension between the HM prior and the cosmological data thus are more favored 
at nearly 3(T than neutrinos with constant mass. These results are shown graphically in Fig[TJ From the left panel 
we see that S is weakly constrained unless we add HM prior while the right panel shows explicitly the modification 
of posterior distribution of neutrino mass today if we allow it to vary with cosmic time. Furthermore, from Figl2]we 
see that the current neutrino mass is correlated with its time variation. This correlation mainly stems from the LSS 



It's true that the bounds on the neutrino mass at current epoch can be relaxed by considering the correlation among neutrino mass 
with dark energy parameters |44|| and with inflationary parameters |45|| . however if the neutrino mass varies with time, its mass limit can 
be much looser. 
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data. We know that the galaxy survey is powerful to weigh neutrinos by detecting the suppression on small scales 
due to the free streaming effect of neutrinos [47l|. The free streaming scale of mass varying neutrinos Xfs-u can be 
roughly estimated as: 



^ 30eV ' 



Mpc 



(10) 



where a^vfl is the scale factor when neutrinos become non-relativistic. We have seen from Egs. lfTU)) and ^ that Xps-u 
is determined by the neutrino mass today and its evolution behavior, which lead to the correlation among m^Q, S and 
TTiiyo * 5. In FigI2]we find that m^o is anti-correlated with S and mj,o * as expected. 





FIG. 2: Contour plots of parameters related to mass varying effect of neutrino, mi is the neutrino mass changing with time, 
say, nil = mo* (5. Different data combinations are distinguished by color. Green: All; Yellow: AU+BAO; Cyan: AU+BAO+HM 
prior. 68% and 95% C.L. contours are illustrated from inside out. 



From Table HI we find that the mean value and the error bars of il^ and erg do not change much if neutrino mass 
varies. At first glance this seems at odds since the neutrino mass limit today has been significantly relaxed by its time 
variation and we know that the neutrino mass today is strongly correlated with matter density and erg as illustrated 
in Fig. ([3]). However since S is anti-corrclatcd with TOj/Qi the aforementioned effect is counteracted, leaves flm and cg 
nearly unchanged. This means mass varying neutrinos can hardly be excluded by data sensitive to Qm and erg, such 
as SN la, CMB and LSS. 

0.4 
0.35 

E 

^ 0.3 
0.25 

0.2 

0.9 
0.85 

0.8 
0.75 

0.7 

12 12 3 

FIG. 3: Contour plots of sum of neutrino mass at present versus Q,m and as. Left panel: Neutrino mass varies as Eq[Tl Right 
panel: Constant neutrino mass. Different data combinations are distinguished by color. Yellow: All; Green: All+BAO; Blue: 
All+BAO+HM prior. 68% and 95% C.L. contours are illustrated from inside out. 
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In Fig.(|4]), we show the evolution of neutrino mass with 1 and 2 a error using all data mentioned in this paper, we 
see that the neutrino mass is best measured at an intermediate redshift rather than now due to the sensitive SN data 
in this redshift range. 

In summary, in this paper, we for the first time study the cosmological implications of the time variation of neutrino 
mass in a model-independent fashion. We find that numerically time variation of neutrino mass can relax the current 



scale factor a 



FIG. 4: Evolution of neutrino mass with respect to the scale factor. From inside out, we show the models with mean parameter 
value (central blue line), 68% C.L.(dark gray) and 95% C.L. (bright gray) fitting with the data combination of AU+BAO+HM 
prior. 



mass limit significantly. This result has some interesting and important implications. It could resolve the tension 
between HM prior, if confirmed, and cosmological data, while it does not spoil the common prediction of Qm and as, 
and it might also be possible to revive the models of warm dark matter which has been shown to be excluded 
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